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ABSTRACT 

We report on sensitive far-infrared observations of ^^CO pure rotational 
transitions in the OMC-1 core of Orion. The lines were observed with the 
Long Wavelength Spectrometer (LWS) in the grating mode on board the 
Infrared Space Observatory (ISO), covering the 43-197 /zm wavelength range. 
The transitions from Jup = 14 up to Jup = 19 have been identified across the 
whole OMC-1 core and lines up to Jup = 43 have been detected towards the 
central region, KL/IRc2. In addition, we have taken high-quality spectra in 
the Fabry-Perot mode of some of the CO lines. In KL/IRc2 the lines are 
satisfactorily accounted for by a three-temperature model describing the plateau 
and ridge emission. The fluxes detected in the high- J transitions {J^p > 34) 
reveal the presence of a very hot and dense gas component (T = 1500 — 2500 K; 
iV(CO)= 2 X lO^^cm"^), probably originating from some of the embedded 
sources previously observed in the H2 near-infrared lines. At all other positions 
in the OMC-1 core, we estimate kinetic temperatures > 80K and as high as 
150 K at some positions around IRc2, from a simple Large- Velocity Gradient 
model. 
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1. Introduction 

The Orion molecular cloud is the prime example of high-mass star forming regions 
and due to its proximity one of the most extensively studied sources. The numerous 
observations performed in the millimeter and submillimeter range (see Blake et al. 1987 
and Genzel & Stutzki 1989 for a review) have shown that the emission arises from three 
main components. The bulk of the molecular gas in the extended ridge (typical densities 
< 10^ cm~^) is heated up to 50-60 K by the UV field of the Trapezium stars. The molecular 
emission near the high-velocity outflow, the plateau^ is distributed in a highly anisotropic 
region of 40" in size with a mean density < 10^ cm~^. The plateau consists of gas heated 
to temperatures of 80-150 K distributed in a low-velocity component centered on IRc2 
and in a bipolar high-velocity component orthogonal to the latter. Finally, the hot core 
(T = 400 — 500 K) of IRc2 is a compact source (10") with densities up to 10^ cm^"^. Its 
high-excitation molecular emission probably arises from material shocked by the interaction 
of the high-velocity outflow with the ambient gas. 

The OMC-1 core has been widely studied in various molecular tracers. Among these, 
CO is of great interest because observations at different transitions delimit regions with 
marked contrasts in their physical characteristics. Most are low- J [J up < 4) observations 
that trace the extended low/moderate density gas, and the bulk of the outflow gas. 
Higher-J transitions (J„p > 6) are sensitive to the denser and higher-temperature gas of 
the plateau or the hot core, however they are much less sensitive to those values from 
the ridge. Multiple CO line observations towards IRc2 carried out with the KAO in the 
submillimeter /far-infrared range (see e.g. Boreiko et al. 1989) suggest evidence for a hot 
gas component at 600-750 K in the high-velocity and (possibly) low-velocity plateau, 
based on the few flux measurements at Jup > 22. However, the large uncertainties in the 
atmospheric transmission yield considerable variations in the parameters of this hot gas 
component. The ISO satellite offered us the possibihty of assessing in a more systematic 
and comprehensive way the temperature distribution in Orion. We present in this Letter 
the flrst large-scale survey of the far-infrared CO lines in the Orion OMC-1 core. 



2. Observations and Results 

Five LWS-grating rasters were performed in the LWSOl mode (range 43-197 /xm) with 
a spectral resolution of 0.3-0.6 iim {R = A/AA ~ 150 — 300). A total of 23 positions (6 scans 
per position) with an angular separation of 90" (~ beam size) were observed around IRc2. 
The central position was included in each raster in order to check the relative calibration of 
the observations. The uncertainties in the absolute cahbration are estimated to be ~ 30%. 
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After pre-reduction (pipeline 7) the data were reprocessed in order to remove non-linear 
effects (Leeks et al. 1999). The size of the LWS beam is actually close to 70" (E. Caux, 
priv. comm.). 

The spectra observed at each position (36 scans for the central one) were averaged 
together in order to minimize the noise level. The statistical errors of the measured lines 
fluxes arc negligible. The transitions at A < 90 //m are effected by additional uncertainties 
in the determination of the continuum level (hence the baseline subtraction), as well 
contamination from neighboring lines. The confusion level is ~ 3 x lO^^^Wcm^^ in the 
averaged spectra. However, we checked that all the features identified in the averaged 
spectra were present in all of the original scans. Figure 1 shows the LWS spectra in the 
range 135-190 //m for the positions where the CO lines could be unambiguously identified. 
Besides CII (158 /xm) and Ol (63 //m and 146 //m) fine-structure atomic hues, the CO 
emission, together with that of H2O, dominates the OMC-1 far-infrared spectrum. Among 
the molecular species detected in Orion at infrared wavelengths, CO appears to be the 
most widespread (a full description of the IRc2 spectrum can be found in Cernicharo et al. 
1999a, 1999b, hereafter C99a, C99b). At all of the positions, we observe the pure rotational 
transitions from J„p=14 to J„p=19, but only the lines up to J„p=17 have a SNR high 
enough to allow confident detection outside IRc2. At large angular distances from the 
central region (> 180"), the CO lines become very weak and the spectra are dominated by 
the atomic fine-structure hues. The CO integrated fluxes at fifteen selected positions are 
given in Table 1. Some of lines are contaminated by nearby features (01 at 145.6 //m, OH 
at 163.2 //m and H2O at different wavelengths), whose contributions have been estimated 
from complementary Fabry-Perot observations (C99a) or from the grating data where the 
line separation was about 0.3 /im. 

Towards the central region, the consecutive pure rotational lines from 7^^=14 up to 
•A(p=43 could be identified. Due to the presence of numerous unresolved components, the 
calibration of the high- J transitions (Jup > 34) is somewhat uncertain. When observing 
the H2O lines towards lRc2 at high spectral resolution with the LWS in the Fabry-Pcrot 
mode (R ~ 7000), several CO transitions, between J„p=14 and 7^^=33, were also detected 
in some detector bands (Fig. 2; for details on the observations, see C99a). The Fabry-Perot 
(hereafter FP) fluxes agree to better than 30% with the grating data, except for the J„p=16 
and the Jup =28 hues where the differences reach 50%. In view of the possible uncertainties 
in the actual pipeline, we have used the grating fluxes for these two lines. Quite remarkably, 
the line intensities are nearly constant from J„p=18 up to J„p=21 indicating that the lines 
are optically thick and that the emitting regions arc similar in size. The line intensity 
merely decreases by a factor of three for the following transitions up to Jup— 28. All of the 
hues have typical widths of Av ~ 60kms~^ (HPFW), i.e. they are partially resolved by the 



-4- 



FP. 

3. Discussion 
3.1. The Ridge 

In order to determine the temperature of the gas traced by the CO hnes detected in 

the OMC-1 core, we have modeled the emission of a gas layer with a column density of 
A^(CO)= 10^^ cm^^ for various densities, n{]i2), and temperatures. The CO fluxes were 
computed using a simple Large- Velocity Cradient approach. For the LWS wavelength 
range the dust opacity is still low enough that the coupling between dust and gas can be 
neglected. The adopted linewidth in the model is At; = lOkms"^ for positions close to the 
center (±90", ±90") and Av = 5kms~^ for the points more distant from IRc2. 

We display in Fig. 3b-d the expected fluxes for the transitions Jup < 25 with n(H2) 
densities in the range 4 — 40 x 10"^ cm~^ for specific temperatures. The temperatures shown 
are those which give the best match to the fluxes measured at three typical positions in the 
core: halfway between the bar and the main core (90", —90"); south of IRc2 towards the 
S6 protostellar source (0", —90"); the extended ridge (0", -180"). For densities in the range 
4 — 40 X 10^ cm~^, typical of those measured in the OMC-1, temperatures larger than 80 K 
are required to account for all of the observed fluxes. Lower temperatures of ~ 60 K fail 
to reproduce the observations by a factor of two or more (see Fig. 3). Around the central 
position, (±90", ±90"), we find evidence for warm gas with T > 120 K and densities of 
the order of ~ 10^ cm~^. The column density adopted in our model appears as an upper 
limit since it was observed towards the IRc2 core (Blake et al. 1987). If we adopt a lower 
column density of 5 x 10^^ cm~^ for positions remote from IRc2, in closer agreement with 
CO observations at miUimeter wavelengths (Bally et al. 1987), reasonable fits are obtained 
by increasing the kinetic temperature by 10-20%, or conversely increasing the density by 
a factor of 2-3. This does not effect our conclusions about the extended warm component 
around the hot core. From Fig. 3, it appears that low-J observations bring very little 
constraints to the temperature determinations in our modeling. The kinetic temperatures 
derived from KAO observations of the Jup = 7 line (~ 90 — 100 K) at 100" resolution by 
Schmid-Burgk et al. (1989) agree reasonably well for the S6 region and close to IRc2. 

On the western side of the ridge, our modeling requires temperatures above 100 K or 
more, and/or high column densities and volume densities to account for the observed fluxes 
(see Table 1). The KAO data suggest moderate temperatures around 50 K, although 
it is difficult to compare both sets of data due to the larger beam and coarse sampling of 
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the former. We note that some contamination in our data from the IRc2 region within a 
possible LWS error beam cannot be excluded and could reduce the derived column densities 
and/or kinetic temperatures. Also the LWS beam profile is known to be rotationally 
asymmetric for an extended source and could account for some of the discrepancies. 

3.2. The Orion-KL/IRc2 Region 

The emission in the far-infrared CO lines from Jup— IS to 33 can be reasonably 
explained by a simple two-temperatures model of the plateau region. First we note that 
a single temperature model cannot reproduce both the flux of the high ( J„p > 28) and 
low-J lines. On the other hand, the range of temperatures reported in the extended ridge 
imply that its contribution becomes negligible with respect to the plateau (see also Fig. 3 
and below) for transitions above J^p = 18. The physical properties of the gas in the 
different regions (in particular the plateau) were chosen as close as possible to those of 
Blake et al. (1987). The low-velocity plateau is modeled as a two-shell region expanding 
at 25kms~^ with a micro-turbulent velocity of lOkms"^. The fine fluxes were calculated 
with the radiative transfer code developed by Gonzalez- Alfonso & Cernicharo (1997) and 
were convolved with the expected instrumental proflle (see Fig. 2) which accounts for the 
FP instrumental response (E. Caux, priv. comm.). We have adopted standard properties 
for the dust (silicate) grains: a typical radius of 0.1 fim and a standard gas-to-dust mass 
ratio of 100. The inner region gives rise to the entire emission of the Jup = 33 and Jup = 28 
lines, to the bulk of the emission of the Jup = 24 line and to a significant contribution of 
the lower- J fines. The colder gas of the plateau also contributes to the emission of the 
low-J fines. Based on our study of the H2O emission in Orion by Cernicharo et al. (C99a), 
we assume that the density is high enough to almost thermafize the CO fines: ~ 10^ cm~^ 
for the inner region and ~ 10^ cm~^ for the external one, close to the estimates derived by 
Blake et al. (1987). We use the standard CO abundance of 1.2 x 10~^. The ratio of the 
Jup = 33 to the Jup = 28 lines sets an upper limit of ~ 500 K to the temperature of the 
inner region: a higher value would imply a lower contrast between the two lines for any 
CO column density. For a lower temperature (~ 350 K), the high column density needed 
to account for the Jup — 33 fine would then overestimate the ^^CO fines, whereas they are 
hardly detected (we adopt an upper limit of 3000 Jy) . We have adopted a temperature 
T — 400 K for the inner plateau region. The resulting column density is A'"(CO)= 10^^ cm~^; 
the outer and inner radii are 3 x 10^^ cm (4.5") and 2.2 x 10^^ cm respectively. For the 
external shell (outer radius = 7 x 10^^ cm), the temperature is T = 300 K and the column 
density A^(CO)= 3.5 x lO^^cm^^. The assumed column densities agree with the values 
quoted in the literature (Blake et al. 1987). 
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The agreement with the hne profiles is very satisfying taking into account the simple 
hypothesis used. The apparent high-velocity emission is fully reproduced by convolving 
with the broad-wing instrumental response (Fig. 2). Down to the sensitivity of these 
observations, we do not detect in this wavelength range any emission from the very 
high-velocity gas around IRc2. Comparison with all the lines observed by the LWS (either 
in FP or grating mode) proves to be satisfying up the J^p = 33 transition. We have added to 
our modehng the contribution of the extended ridge, as a layer at T = 80K (see Sect. 3.3) 
with a column density of A'^(CO)= 4 x 10^^ cm~^. The resulting fit and the contributions 
of the different regions are shown in Fig. 3. As a whole, the emission observed towards the 
IRc2 core can be satisfactorily accounted for by two regions of the plateau gas at different 
temperatures: ~300 K and ~400 K, and from the ridge contributing significantly up to the 
Jup = 16 line. Note that since the central part of the line profiles is only partially resolved 
by the FP, we cannot exclude that a part of the warm component modeled actually arises 
from the hot core (where At' ^ 5kms~^). Therefore, the emission from the T = 400 K 
component should be considered as the combined emission of the plateau and the hot core. 
The contribution of the ridge seems to be underestimated in our model. However, as for the 
adjacent positions, the contribution from an error beam could effect the observed fluxes. 

3.3. The High-J CO lines 

The higher- J transitions (J„p > 34) observed with the grating mode are largely 
contaminated by other emission lines. Using our FP data we have estimated the contribution 
of the strongest adjacent fines (OH, H2O and NH3) to the CO fluxes. In view of the high 
density of spectral lines for A < 90 //m, additional contamination by other weaker lines 
cannot be discarded (C99a,b). We note, however, that the lines are well above the confusion 
level and that most of the observed flux must arise from CO itself. 

After correction, all the CO lines observed either in FP or grating mode lay above the 
model of the plateau and the ridge (Fig. 3a), showing evidence for a hotter gas component. 
In the absence of other information on the spatial location and the kinematics of this 
gas, we speculate that such emission arises from the shocked regions of interaction of the 
Orion outflows with the ambient medium, detected in the 2.12/im lines. All sources of 
strong H2 emission (BN, KL, PKl, PK2; Beckwith et al. 1983) are located within our 
70" beam. We modeled the high-J CO emission using a Large- Velocity Gradient (LVG) 
approach with parameters representative of the H2 Peak 1 region: a size of 10", a density 
n(H2) ~ 3 X lO^cm"^ for the shocked emitting region, and a finewidth Av — 30kms~'^. 
The fluxes are roughly accounted for by a gas layer at T ~ 1500 — 2000 K and a column 
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density A'"(CO)~ 6 x lO^^cm"^, implying a shell thickness of 1.6 x 10^^ cm. These values 
are merely indicative due to the lack of spectral and spatial information, but are consistent 
with those predicted by shock models (e.g., Flower & Pineau des Forets, 1999). 
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Tabic 1 : CO line fluxes measured in the 
extended ridge of OMC-1. 

The positions for which the line fluxes could be contaminated by 
emission from IRc2 within the LWS error beam are marked with an 
asterisk (see text). 
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Fig. 1. — LWS grating spectra from the observed positions around IRc2 between 135- 
197 /im. Lower marks indicate the ^^CO rotational transitions from J^p — 14 to J„p = 19. 
At the central position ^'^CO transitions are marked by the arrows. Some spectra have been 
multiplied by a scale factor indicated at the top right of each panel. 

Fig. 2. — A montage of select CO transitions observed towards IRc2 with the LWS in the 
Fabry- Perot mode (thick hne). We have indicated (thin line) the modeled emission of the 
plateau (see text). The instrumental response of the Fabry- Perot at 145 //m is indicated on 
the — 18 spectra (dotted line). 

Fig. 3. — a) Integrated flux of the CO lines detected towards the IRc2 region. The modeled 
contributions from the plateau (P, thin line) and the ridge (R, dotted thin) are shown. 
A hot gas component (HG) was introduced (dashed line) to reproduce the emission at 
Jup > 34. The arrows mark the fluxes of the observed CO lines possibly overestimated 
due to contamination by weak adjacent lines. We have superposed (empty squares) the 
values estimated for the ISO beam size from previous work: Jup= 4 (Schultz et al., 1992); 
6 (Graf ct al., 1990); 7 (Schmid-Burgk ct al. 1989; Howe ct al. 1993); 21 (Boreiko & Betz, 
1989); 22-34 (Watson et al 1985). b,c,d) Expected fluxes for the GO lines up to Jup = 25 
observed in the extended ridge with the ISO LWS beam, for kinetic temperatures comprised 
between 60-120 K, H2 densities 4 - 40 x 10^ cm-=^ and a column density A^(CO)= 10^^ cm-^. 
We show the "best temperature" fit, for each pointing, which accounts for the measured 
fluxes (Table 1). The arrows have the same meaning as above. 
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